A test of the (7— ology model for one-dimensional interacting Fermi systems 
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Bosonization predicts that the specific heat, C{T), of a one- dimensional interacting Fermi system 
is a sum of the specific heats of free collective charge and spin excitations, plus the term with the 
running backscattering amplitude which fiows to zero logarithmically with decreasing T. We verify 
whether this result is reproduced in the p— ology model. Of specific interest are the anomalous 
terms in C{T) that depend on the bare backscattering amplitude. We show that these terms can be 
incorporated into a renormalized spin velocity. We do this by proving the equivalence of the results 
for C{T) obtained within the g-ology model and by bosonization with velocities obtained by the 
numerical solution of the Bethe-ansatz equations for the Hubbard model. 

PACS numbers: 



Onc-dimensional interacting fermionic systems are be- 
lieved to be described by an effective low-energy theory: 
the g-ology model flSd. This model involves a small 
number of interaction vertices describing small momen- 
tum scattering of fermions near the same and opposite 
Fermi points (34 and (72, correspondingly) and 2kF scat- 
tering gi [for fermions on a lattice, there is also an Umk- 
lapp vertex (33)]. The effective vertices are, in principle, 
obtained by integrating out high-energy fermions in mi- 
croscopic models, e.g., the Hubbard model. To first order 
in [/, gi — U/{2'KVp); beyond first order, all gi are differ- 
ent, and each of them is represented by a series in U. 

A powerful way to treat the g-ology model is bosoniza- 
tion, which transforms interacting fermions into the col- 
lective bosonic excitations in the charge and spin chan- 
nels 0, Q . For the case of a repulsive interaction between 
original fermions, considered in this paper, the bosoniza- 
tion shows that the charge sector is a free Gaussian the- 
ory, while the spin sector becomes asymptotically free in 
the low-energy limit, when the coupling of the marginally 
irrelevant process (2fci? scattering of fermions of oppo- 
site spins) flows to zero d, 0, H, Q As a consequence, 
bosonization predicts that at the lowest temperatures the 
specific heat of interacting ID fermions is the same as the 
specific heat of two systems of acoustic f D phonons, i.e., 

CiT) = ^ (- + -) , (1) 

where Vp and Va- are the charge and spin velocities, cor- 
respondingly. Eq. II]) was verified by a weak-coupling 
renormalization group (RG) treatment of the original 
fermionic model, in which all non-logarithmic corrections 
to the couplings were neglected [ITSBIi^l- However, 
it has never been proven explicitly that the bosonization 
result in Eq. ([T]) is valid to all orders in the interaction. 

Recent results call for a further study of the validity 
of Eq. ([1]) beyond the weak-coupling limit. In particular. 



two of us have obtained [11] the specific heat of interact- 
ing ID fermions to order gf and found that the effect of 
2kp scattering in the spin channel is more involved than it 
had been previously thought - in addition to terms that 
depend on the running coupling gi(T), the free energy 
also contains terms that depend on the bare coupling gi 
(see Eq. ^ below). 

^From the field-theoretical point of view, these terms 
should regarded as anomalies, i.e., they can be viewed 
equivalently either as low- or high-energy contribu- 
tions. These anomalous terms cannot be simply absorbed 
into the Gaussian part of the bosonic Hamiltionian, and 
it is not a priori clear whether such terms can be incor- 
porated into the renormalized charge and spin velocities. 

The goal of this paper is to prove that the answer to the 
question formulated above is affirmative, at least within 
perturbation theory. We show this by comparing the spe- 
cific heat obtained from second-order perturbation the- 
ory in the g— ology model, calculated in Ref. Il', with 
Eq. ([1]) , where the velocities are obtained from the Bethe 
ansatz solution of the Hubbard model [4]. 

Our starting point is the Hamiltonian of the ID Hub- 
bard model 

where ek — ~2tcosk and /i = —2t cos hp is the chemical 
potential (the lattice spacing is set to unity). For U — 
the Fermi velocity is related to the Fermi momentum hp 
by Vp — 2t sin kp. 

Away from half-filling, umklapp scattering requires 
collisions of more than two fermions and is therefore ne- 
glected in our treatment. The spin and charge veloc- 
ities can be obtained by solving certain linear integral 
equations obtained from the Bethe ansatz solution of the 
Hubbard model [4] . At strong coupling, the spin velocity 
is vs = (7r/2)(4<V[/)(l - sin2kp/2kF) [13]. The weak 
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TABLE I: Coefficients c and c for different densities in tfie ID 
Hubbard model. Tfie density n = 2kF/'T; n — 1 corresponds 
to lialf-filling. 



coupling limit has been studied in [l3|. It is possible to 
derive a small U expansion for the spin and charge veloc- 
ities analytically by solving the integral equations drived 
in Ref. by Wiener-Hopf methods. As we are inter- 
ested in the 0{U'^) terms of the expansions, the resulting 
calculations are somewhat involved. We therefore have 
solved the integral equations numerically for small U and 
found that the results are well fit by the series 




(3a) 
(3b) 



Based on the analytic result for the spin velocity at half- 
filling [3] and the fact that the only marginally irrelevant 
operator (the interaction of spin currents) is the same at 
and below half-filling we do not expect logarithmic terms 



to appear in these expansions. We have determined the 
coefficients c and c' for a number of different densities. 
The results are shown in Table U As we will see below, we 
will only need the difference c' — c to verify the validity 
of Eq. IH) . In all cases we find that within the numerical 
accuracy of our computation 



(4) 



The deviation from 1 is less than one percent in all cases 
The spin and charge velocities can alternatively be calcu- 



lated in perturbation theory for the 5— ology model [11 
141 . An explicit calculation of C{T) within this model 
is somewhat involved as the low-energy g— ology model 
contains two momentum cutoffs, A^ and A;,, constrain- 
ing the integration over the fermionic dispersion and over 
the momentum transfers near 2k p, respectively. [The 
g-ology model is only valid when Af > Ab, i.e., when 
the interaction vanishes at the cutoff set by the disper- 
sion.] Some terms in C{T) are cutoff-independent while 
some depend logarithmically on the ratio Af/Ab. Fortu- 
nately, at least to second order in g, all cutoff-dependent 
renormalizations can be absorbed into the renormalized 
backscattering amplitude gi = gi — log (Af/Ab), so 
that the specific heat is expressed in terms of (74, g2, and 
gi without any explicit dependence on the cutoffs [ll| . 
To second order in g, C{T) is given by 



C{T) = 



2-kT 
3vf 



1 + (51 - 54) + (51 - gif +gl + 
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0(9') 



(5) 



where all vertices are measured in the units of 2ttvf- 

To compare Eq. ([5]) with Eqs. ^ where Vp and Uo- 
given by Eqs. pa|3bp . we first note that the g-ology 
model can be bosonized by expressing the operators of 
right- and left-moving fermions, Ra and {a —"l, |) as 

Ra{x),La(x) = ^_ exp [±i {(t)a {x) T Oa{x))] , (6) 
v27ro 

where & is a short-distance cutoff related to the fermionic 
momentum cutoff (Af) of the t/— ology model. Under 
bosonization, the terms in the fermionic Hamiltonian pa- 



rameterized by the couplings 54 and 52 are mapped onto 
the free, Gaussian part of the bosonized Hamiltonian. 
The 2k F term, parameterized by gi, leads to non-linear, 
cosine terms in the bosonic Hamiltonian, which give rise 
to interactions in the spin channel. 

To first order in gi, backscattering just renormalizes 
the prefactors in the Gaussian part of the bosonized 
Hamiltonian [3, so that the 17— ology model can 
be reduced to a gas of free acoustic bosons with Hq ~ 
H^G^ +H^J\ where 



h: 



(p) 



^ I dx{l + 254 + 252 - 2gi) [d^c^pf + (1 + 254 - 252) {djpf , 



^ = ^ / cfx (1 - 251) [d^cj^^f -t- {d^e^f , 



(7) 
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FIG. 1: One-loop diagrams for the interaction vertices. In 
this paper, we need only renormalizations coming from the 
high-energy scales. The high-energy renormalizations of g4 
and gi come only from the Cooper diagram d). 



and the charge and spin bosons are defined as 4>p,cr — 
(0T ± 0i) /V2 and 0p,, = (0T ± ^i) /V2. 

If this were the only effect of backscattering, the spe- 
cific heat would be given by Eq. ([T]) with the effective 
spin and charge velocities Va and Vp read off from Eq. 
0: 



il^vj,{il + 2g4~gif ~ (232-91 f) 



(8) 



Re-expressing Eq.([5]) in terms of Va and Wp, we obtain 



C{T) 



ttT 
~3" 



TiT 



- + -) + 3^3f + (9) 



We see that this expression differs from Eq. ([T|). This is 
not surprising because the effect of backscattering can- 
not be simply absorbed into the Gaussian part of the 
bosonized Hamiltonian, beyond the first order in gi. 

The issue therefore is whether the extra gl term in Eq. 
^ can be absorbed into renormalization of the veloci- 
ties v ^ so that the specific heat is still given by Eq. 
([1]) with the renormalized velocities Vp and Va . A simple 
extension of the previous analysis to a non-SU(2) sym- 
metric case shows that, beyond the first order, 2k p scat- 
tering contributes only to the spin part of the bosonized 
Hamiltonian. The real issue then is renormalization of 
the spin velocity Ug- — > Wo-. The charge velocity given by 
Eq.Q must be the same as the exact one, i.e., Vp = Vp. 

A straightforward way to check this is to compare Eqs. 
dl]) and dHI with Eq. pa|3bp . Quite generally, one can 
write Va- — Va- + avpgl, where a is a dimensionless con- 
stant. Only if a — —1, the extra gf term in Eq.® can 
be absorbed into Va, i.e., Eq. ^ can be cast into Eq. 
((H). Using Eq. dS]), we obtain 

Vp ^vf[i + 254 - 51 - ^(252 - 5i)' j , (10a) 



Va=VF [I - gi- 



1 - 2a 



-9l 



(10b) 



These two expressions should be the same as Eqs. pa|3bp . 

To compare Eqs. Pal3bp and (|10a|10bp . we need to 
evaluate renormalizations of 54 and 51 to second order 



in U and to select the contributions which comes from 
high energies. By construction, such contributions are 
absorbed into the bare couplings of the 5— ology model. 

There are four second-order diagrams for the 54 ampli- 
tude to order U'^ (see Fig. 1). For a constant U, diagrams 
a) and b) cancell each other. Diagram c) contains the po- 
larization bubble n(0) — (l/27r)^/ dkduj/{iuj — ek+ fJ-)^- 
Renormalizations due to 11(0) are within the low-energy 
theory, as one can evaluate 11(0) in such a way that the 
result n(0) = — l/(7rt>_F) is determined entirely by the 
states near the Fermi energy [ll|. The remaining dia- 
gram d) describes renormalization in the Cooper channel. 
Up to a prefactor, it is given by 



dkduj- 



1 



The momentum integration in Eq. (jlip is not confined 
to the Fermi surface, i.e., this diagram does contribute 
to high-energy renormalization of 54. In the Hubbard 
model, the momentum integration is limited from above 
by the Brillouin zone. The lower limit Af can be safely 
set to zero as the integral is infrared- finite, i.e., J dk = 
dk. Integrating over lu and then over k, we obtain 



U 



94. 



2'kvf 



1 



U 



2nvp J 



(12) 



Note that the only depence on the density is through 
Vp — 2tsmkF. We verified that a model of fermions in a 
continuum with dispersion k"^ /2m gives the same result. 

Renormalization of 51 is more involved because 51 is 
a running coupling, and the second-order result depends 
logarithmically on the upper cutoff of the low-energy the- 
ory, which is the lower limit of the integration for "high- 
energy" renormalization. The diagrams are the same as 
in Fig. 1. As for 54, diagrams a) and b) cancel each other 
while diagram c) contains H(0). Therefore, only Cooper 
diagram d) contributes to the "high-energy" renormal- 
ization. Evaluating Cooper diagram for backscattering, 
we obtain 



51 = {UI2t:vf) 1 



U 



2'KVf 



where 



Up „ 
L = —V 
27r 



dk 



duj- 



1 



(13) 



(14) 



and The symbol V J here implies that momentum inte- 
gral does not include the regions near the Fermi points 
of width 2Af. Integrating over frequency and then over 
momentum, we obtain 

2 sin kp 



L = 21og- 



Af 



(15) 



Substituting Eq. into Eq. p^ . we then obtain 

U / U 2 sin kp 
91 = 7, 1 log — I 

ZTTVp \ TTVp Af 



(16) 
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In addition, gi is renormalized within the low-energy g- 
ology model. This renormalization depends logarithmi- 
cally on the ratio of the bosonic and fermionic cutoff 
of the g-ology model: gi = gi - log (Af /Ab) [ll|. 
Adding up this result with Eq. p?]) we find that the com- 
bination of the high-energy and low-energy renormaliza- 
tions just replaces Af by Ab under the logarithm, i.e., 

C/ / C/ , 2smkF\ 
91 = o 1 log —A ■ 

Alternatively, one can obtain the full renormalization 
of gi by excluding the regions of width Ab near zL2kF 
from the integration over the momentum transfer k. This 
gives the same result as in Eq. (IT71) . 

The value of Ab is unknown: as we said earlier, the 
ology model assumes that the low-energy properties 
of the original system of ID fermions with a short-range 
interaction are the same as in the model where interac- 
tions are artificially restricted to narrow regions of mo- 
mentum transfers either near zero (for g^ and 32) or 2k p 
(for gi). We can only realistically expect that vp^h is 
substantially smaller than a half of the fermionic band- 
width W/2 = 2t. Still, we have two pairs of equations to 
compare [Eqs. ()3a|3bp and ()10a|10bp ] and two unknown 
parameters: Ab and a. Solving for the unknowns, we 
obtain 

Ab = 2sinfcF ''^p(^^ ~ I) ' 
a = c'-c-2. (18) 

By virtue of Eq. we conclude that a — ~1, which is 
precisely the value of a one needs to cast Eq. ([9]) into 
Eq. ID). This, we believe, is a "numerical proof" of the 
statement that at very low temperatures the specific heat 
of ID interacting fermions is the same as two system of 
acoustic phonons with certain spin and charge velocities. 

We also find that for most of densities vp-^b/i'^t) — 
Af, sin kp is smaller than one, as it should be, otherwise 
the g— ology model cannot be justified. In particular, at 
quarter-filling, Abs'mkp ~ 0.44. Near half-filling, how- 
ever, Absinkp becomes larger than one, which questions 
the validity of the 5— ology model. Note also that in the 
opposite limit of small density. A;, goes to zero as it in- 
deed should as at vanishing kp the linearized dispersion 
no longer holds. 

To conclude, in this paper we obtained expressions for 
spin and charge velocities for interacting ID fermions in 
terms of the couplings of the g-ology model. These can be 
related to the microscopic parameters of the ID Hubbard 
model via a comparison to weak coupling expansions of 
the velocities obtained from the Bethe ansatz solution. 



Using these results, we have shown that all terms in the 
specific heat in the ology model that do not flow under 
RG are absorbed into the specific heat of two free gases 
of massless bosons. As a result, the full specific heat of 
a system of interacting fermions in ID is a sum of the 
specific heats of two free massless Bose gases of charge 
and spin excitations, and the true interaction term, which 
contains the running backscattering amplitude and loga- 
rithmically flows to zero with decreasing T. 
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